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O ' Abstract 

(N : 

^ ^ We consider two fractional versions of a family of nonnegative integer valued processes. We 

■ prove that their probability mass functions solve fractional Kolmogorov forward equations, and 

we show the overdispersion of these processes. As particular examples in this family, we can 
define fractional versions of some processes in the literature as the Polya-Aeppli, the Poisson 
I Inverse Gaussian and the Negative Binomial. We also define and study some more general frac- 

tional versions with two fractional parameters. 
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1 Introduction 

In this paper we consider a large family of nonnegative integer valued processes {M{t) : t > 0} 
defined by 

^ ■ ■■= E (1) 

. where {X„ : n > 1} is a sequence of i.i.d. positive and integer valued random variables, independent 

of a (non- fractional) Poisson process {N\{t) : t > 0} with intensity A. Throughout this paper 
we also deal with a mixed Poisson representation of the process in ([1]); more precisely we mean 
{Ni{S{t)) : t > 0}, where {S{t) : t > 0} is a subordinator, independent of a (non-fractional) 
Poisson process {Ni{t) : t > 0} with intensity 1, such that Ni{S{t)) is distributed as M{t) for each 
fixed t > 0. 

. Our aim is to present and analyze two fractional versions of the process {M(t) : t > 0} (see 

(FVl) and (FV2) at the beginning of Section [3|). This allows to define fractional generalizations 
of some processes in the literature, which include the Polya-Aeppli, the Poisson Inverse Gaussian 
and the Negative Binomial. These processes are commonly used when the empirical count data 
exhibit overdispersion, i.e. when the sample variance is larger than the sample mean; moreover it 
is known that mixed Poisson processes provide simple counting models with overdispersion. We 
shall see that the fractional versions have the same feature. 

The two fractional versions are obtained by considering independent random time-changes of 
{M{t) : t > 0} in terms of a stable subordinator A"' or its inverse (for a G (0, 1)). This approach 
is inspired by the recent increasing interest on random time-changed and subordinated processes; 
see, among the others, [31], [22], [7], [18] and [16]. These processes are widely studied and applied, 
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mainly because their finite dimensional distributions display power law decay and thus heavy tails, 
either in the case of renewals and birth processes (see [20j and f25] ) or in that of diffusions (see , 
[29j and [24j). We recall that fractional versions of the compound Poisson process (as the ones in 
this paper) have been studied in [32] and [5], in the case of continuous distributed summands. 

In this paper we provide a natural fractional extension of some results in the literature. In 
particular Proposition 13.11 below shows that the probability mass functions of the two fractional 
versions of the process in ([T]) solve suitable Kolmogorov equations, where the classical derivatives 
are replaced by the fractional derivative in the Caputo sense (for the first version), and by the right 
sided fractional Riemann-Liouville derivatives on M-|_ (for the second version). 

We conclude with the outline of the paper. We start with some preliminaries in Section [2j In 
Section[3]we illustrate some properties of the fractional versions of the process in ([1]). Some examples 
are presented in Section HI Finally Section [5] is devoted to more general fractional versions of some 
processes in Section [J] (the Polya Aeppli and the Poisson Inverse Gaussian) with two fractional 
parameters. 



2 Preliminaries 



We start with some preliminaries on fractional calculus, i.e. we give the definitions of two fractional 
derivatives for real functions defined on [0,oo). 

1. If € (0, 1), the Caputo derivative of order v (see e.g. (2.4.17) in |llj with a = 0) is defined 

by 

Remark. If we set 1^ = 1, cDq_^ ^ coincides with the classical derivative ^; see e.g. Theorem 
2.1(b) in m]. 

2. If € (l,oo), the right sided fractional Riemann-Liouville derivative on M_|_ of order v (see 
e.g. (2.2.4) in [11]) is defined by 

1 / d\'^ r°° f(s) 
-°VW-f(;;7ri;)(-5) / (f. alU > O), whe.s ™ := h + l. 

Remark. If we set u = 1, rlD^ i coincides with — ^ (i.e. the opposite of the classical 
derivative); see e.g. (2.2.5) in [TTj . 

Hereafter, for simplicity, we always write 

^ := cDl^^t (for u e (0, 1)) and ^ := rlD-_^^ (for v G (1, oo)). 

Throughout the paper we often use the symbol Z{-) to mean a process {Z[t) : t > 0}. In view 
of what follows we recall some preliminaries on the stable subordinator .A"(-) of order a G (0, 1), 
and its inverse ^°'{-)- More precisely let A^{-) be the Levy process (starting at the origin) such 
that, for each fixed t > 0, we have 

Er-W] = I -p(-(-^)"^) ii^i (2) 

^ ^ [ oo if 6* > 0; ^ ' 

thus, by referring to [30j . the random variable A°^{t) has stable distribution of index a and param- 
eters /i = 0, 6* = 1 and a = (t cos(^))^/°. Furthermore £"(•) is defined by 

£"(t) := inf{z : A'^{z) > t}; 
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then we have P{L"{t) < z) = > t) for all z,t > and 

]E[e-e£"(t)] = Ea^i{-ee) for all > 0. (3) 

Now we recall the definition of the two following fractional Poisson processes (see e.g. [21] for 
u € (0, 1) and [26j for the case u € (1, oo)). In both cases we consider a random time-change of a 
non- fractional Poisson process N\{-). 

1. For u G (0, 1), Nl{-) := NxiL^i-)) where Nx{-) and are independent. 

2. For u G (l,oo), N^{-) := NxiA^^^i-)) where Nx{-) and are independent. 

We can give some further details and we need to introduce some notation. For all integer r > 
and for all 7 G M, the rising factorial, also called Pochhammer symbol, is defined by 



the falling factorial is defined by 

(7). ■■= 

and we also consider the notation 



7(7 + 1) • • • (7 + r - 1) if r > 1 
1 if r = 0, 



7(7 - 1) • • • (7 - - 1)) if r > 1 
1 if r = 0, 

(7)r _ 7(7-l)---(7-r+l) 



r J I 1 if r = 0. 

1. If G (0,1), it is known (see [21j) that N^{t) is distributed as Yln>i^{Ti-\ \-Tn<t}j where 

{Tn : n > 1} are i.i.d. random variables with Mittag-Leffler distribution, i.e. with continuous 
density 

f{t) = Xe-'E,,,{-Xni(^o,oo){t) 
(see e.g. |3j and [H]), where Ea^p{x) := Z^r>o r(ar'+/3) Mittag-Leffler function (see 

e.g. [28], page 17). Then, if we consider the generalized Mittag-Leffler function £'2^(x) := 

T,r>o rirlar+fi) ' '^^^^ (^^^ formula (2.5) in [4J) 

P{N^{t) = k) = {Xf'fEl+l^^^-Xt'') for all integer > 0. (4) 

2. If z/ G (1, 00), we have 

P(N^(t) = k) = ^—^ V ^ r^(r/iy)k for aU integer A; > 0. (5) 

ki ^-^ r 

r=0 

We recall that this process is presented in [26] by referring to the fractional difference operator 
(1 — B)", for a G (0, 1], and it is called space-fractional Poisson process. More precisely the 
probability mass function in ([5]), with = 1/a, satisfies the following equations: 



|P(jyV«(^) = /c) = -A°(l - B)°'P{Nl^"{t) = k) for ah integer k>0 
P(iV|/"(0) = 0) = 1, P(iV|/"(0) = k) = for all integer k > I, 



where B is the so called backward shift operator defined by Bf{k) = f{k—l) and B^ ^Bf{k) 
f{k — r); thus in particular we have 



00 



[i-Brf{k) = ^{-iy[ ]fik-j). 

j=0 



a 



We remark that formula ([5]) here coincides with formula (1.2) in [26] . though it is written in 
a slightly different way. 
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3 The fractional versions: compound and mixed representations 

We start with the compound representation of the fractional versions of the process M(-) in ([1]), 
and we refer to the processes defined in Section [2j 

(FVl): M-(t) := M{L'^{t)) = ^k, for u e (0, 1); 

(FV2): M-(t) := M{A^/''{t)) = Y.T=? ^fc, for ^ G (l,oo). 

In view of what fohows we consider the following notation for the probability mass function and 
the probability generating function of the random variables {Xn : n > 1}: 



qk ■= P{Xi = k) (for all integer k > 1); gqiu) '■= u'^qk- 

k=l 

We also consider the notation 

qr:=P{Xi + ...+Xn = k), (7) 

and therefore we have f;^" = for all integer k < n. 

Now we introduce the two following probability mass functions: 

1. for v G (0, 1), we have 

v'^U) - PiM^it) -k)-{ ^(5^') = °) ^ = ° f8l 

p,[t) .- P{M, (t) - fc) - I ql-P{N^,it) =n) if > 1; 

2. for z/ € (1, oo), we have 

I En=i % P\N\{i) = if A; > 1. 

The aim of this section is to prove that these probability mass functions satisfy suitable versions 
of Kolmogorov equations with fractional derivatives. 

Proposition 3.1. If v ^ (0,1), the probability mass function in ([8]) satisfies the equations 
( dy„ 



I m^Plit) = -Xplit) + AEti qiPl-i{t) for all integer k>l, 

with the initial conditions Pq{0) = 1 and p^(0) = for all integer k > 1. 
If v ^ (l,oo), the probability mass function in ([9]) satisfies the equations 

I ^Pkit) = ^Plit) - AEti QrPUit) for all integer k > 1, 



(10) 



(11) 



with the initial conditions Pq{0) = 1 and j5^(0) = for all integer k > 1. 

Remarks. For v = 1, both formulas (|10p and (|lip reduce to the well-known Kolmogorov forward 
equations (for instance one can specialize eq. (2.3) in Chapter 14 in UO^ )- For k = 0, the equations 
in and ([TT]) coincide with the well-known fractional relaxation equations (see, for example, 
l^). For k > 0, they can be seen as a discrete version of the fractional master equation (see eq. 
(5.14) in f21\l for the case where the common distribution of the jumps have continuous density; 
see also JB^). 
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Proof. In both cases {v £ (0, 1) and v G (1, oo)) the initial conditions trivially hold. 

Case V € (0, 1). We have to check the equations in (jlOp . It is known (see e.g. Theorem 2.1 in [4j) 

that 

f ^P(iVj:(t) = 0) = -\P{N'^{i) = 0) 

I J^(^a(*) = k) = -XP{N^{t) = k) + \P{Nl{t) = k-l) for all integer k>l, 

with the initial conditions 

P{Nl{<d) = 0) = 1 and P{Nl{{)) = /c) = for aU integer k>l. 

The initial conditions for (P(A^j^(0) = h))^^^ meet the ones for (p^(t))fe>o in the statement of the 
proposition. For A; = we have ^Po(i) = —XPo{t) because Po(t) = P{N^{t) = 0). For A; > 1 we 
have several steps. Firstly, since 

h 

pl{t) = ql^'PiN^it) = n) for ah integer h>k, 

n=l 

by ([El) we get 

k 

^pUt) = J2 an-mm) = n) + XP{Nm = n-l)] 

n=l 



k k / k \ 

n=l n=l \i=l / 

k k 

Xplii) + XY,<l^iZ itr'^'Pim) = n - 1); 



n - 1 



( 

i=l n=l 

moreover, since qf = l{h=o} and Qq'' = l{h=o}, we get 

E E %t-'^pimt) = n - 1) = E 9. E %t-'^p{Nm = n-i) 

1=1 ra=l i=l n=l 

+ftE«o^""'^^™*) = ^-i) 

71=1 

k—1 k 

= E «^ E llt^'^Pi^'xit) = n - 1) + q,PiN^,{t) = 0) 

i=l n=2 

k-1 k~l 

= E E 1k-^PiNxit) = j) + iMt) 

i=l i=l 

fe-1 k 

= E liPk-iit) + QkPoit) = E QiPk-iit)^ 
i=l 1=1 

and this completes the proof. 

Case € (1, oo). We have to check the equations in (jlip and we follow the same lines of the previous 
case. The main difference concerns only the initial step, i.e. the analogue of (jl2p presented above. 
Thus we only give some details on how to prove that we have the following equations 

r ^p(^N^^^t) = 0) = XP{N^^{t) = 0) 

\ ^P{Nl{t) = k) = XP{N^{t) = k)- XP{Nl{t) = k-l) for all integer A; > 1, 
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with the initial conditions 

P{N^{0) = 0) = 1 and P{Nxi^) = k) = for all integer k>l. 
Thus, since N^{-) = NxiA^^" {■)), we have 

.\k 

-e^^^ f,y{x,t)dx for all integer k >0, 



PiN^t) = k) 



kl 



where f^{x,t) := fjii/u{t)i^) density of the random variable A^/'^{t). Then we get (the first 

equality holds by eq. (5.17)-(5.18) in [5] with 7 = ■^) 



"'RL 



(Ax) 
k\ 



-Xx 



fu{x,t) 



x=0 



{-\P{Nx{x) = k) + XP{Nx{x) = k-l)) Uix,t)dx 
=XP{Nm = k)- \P{Nl{t) = k-l), 
and the proof is complete by taking into account P{N^{t) = —1) = for the case A; = 0. □ 

Remark. If € (1, oo), we can also obtain the alternative equation 
d 



dt' 



-Mt) 



-X^'''Y,{-^y{^^''^ ^if^^^Pi^xit) = I) for aU integer k>l, 

j=0 ^ ^ 1=0 



by taking the classical derivative with respect to t in j5^(i) = Yln=i iV^^i^xi'^) — by ([6]). 

It is easy to check that, for = 1, we meet the classical Kolmogorov equation. 

Now let A'^i(-) be a non- fractional Poisson process with intensity 1, independent of £'^(-) (for 
u G (0, 1)) and of yi^/'^(-) (for g (l,oo)). Moreover let S{-) be a subordinator, independent of all 
the other processes such that the following conditions hold. 

1. If z^ G (0, 1), Wit) is distributed as A^i(5(£^(t))); in what follows we set S^{-) := S{L^{-)), 
and we have 



k\ 



E 



k\ 



for all integer /c > 0. 



(13) 



2. If z/ G (l,oo), M^(t) is distributed as NiiSiA^/"" {t)))] in what follows we set 5^(- 
5(yiV'^(.)), and we have 



E 



(5(yiv-(t)))^ ,^(^:/.(,^^ 

A:! 



for all integer A: > 0. (14) 



In such a case we talk about of mixed representation of the fractional versions of the process M(-) 
in ([1]). Furthermore, in view of what follows, it is useful to introduce the function ks defined by 
Ks{6) := logE[e^'^^^^] for all G M (note that ks{-) is a nondecreasing function and ks{0) = 0). 
This function has a strict relationship with the probability generating function gq presented above; 
actually, by considering standard computations on independent random time-changes for Levy 
processes (see e.g. [6]) one can easily check that 



E 



JSit) 



,Xt{gq(l+e)-l) 



(15) 
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for all 6 such that 1 + 9 belongs to the domain of convergence of gq. 

We conclude with the following further formulas: if € (0, 1), by ([3]) we have 



E 



if G (1, oo), by ([2]) we have 



E 



E^^i{Ks{-dy) if 0>O; 



exp(-(-K5(0))^/^t) if < 
oo if 6* > 0. 



(16) 



(17) 



On the concept of overdispersion. It is well-known that a real valued random variable Y is 
said to be overdispersed if Var[y] — E[y] > 0; similarly, for a process Y{-)^ we have overdispersion 
if all the random variables {Y{t) : t > 0} are overdispersed. Typically the compound Poisson 
process M(-) in ([1]) exhibits overdispersion when we exclude the trivial case where the random 
jumps {Xn : n > 1} are all equal to 1; actually, in such a case, we have M(-) = Nx{-), i.e. M(-) is 
a non-fractional Poisson process, and we have Var[M(i)] — E[M(t)] = for all t > 0. 

Here we want to study the same feature for the process M'^(-) in (FVl); we do not deal 
with M'^(-) in (FV2) because each random variable M^{t) has infinite mean (actually the moment 
generating function oiA^/'^{t) is not finite in a neighborhood of the origin; see formula ([2]) presented 
above). It is known (see e.g. formulas (2.7) and (2.8) in [3j) that 



mm] = 

Var[iVj:(t)] 



_ Xt" 



+ {Xt^fZ{u), where Z(i/) := i 



Thus, since we have Zlu) > for all E (0, 1) and Z(l) 
process N^{-) (for u € (0, 1)) exhibits overdispersion. 



(see Figure 1), the fractional Poisson 




Figure 1: Plot of Z{i^) : = 



1 / 1 



versus v. 



Furthermore, it is known (see e.g. Appendix B.4 in [12]) that 

E[M'^(t)] = E[7Vj;(t)]E[Xi] 

Var[Af^(t)] = E[iVj;(t)]Var[Xi] + Var[iVj;(t)]E2[Xi]; 
then, with some computations, we get 

Var[M''(t)] -E[M''(t)] =E[Nl{t)]{\ax[Xi] -¥.[Xi]) +\iXT[Nl{t)]¥?[Xi 



mx't] - E[Xi]) + {Xt''fZ{ij)E'[Xi] 
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and — E[Xi] > because P{Xi > 1) = 1. In conclusion the compound process M^(-) in 

(FVl) exhibits overdispersion. 



4 Examples for the process M( ) 

In this section we present some examples for the process M(-) in ([1]) which allow to define the 
corresponding fractional versions M'^(-) in (FVl) and M^{-) in (FV2). More precisely we consider 
suitable choices of {qk)k>i and A; moreover, in all cases, the probability mass function {qk)k>i is 
a zero truncated negative binomial distribution as in Example B.3.1.5 in Appendix B.3.1 in [12] . 
where /3 in that reference stands for for some a € (0, 1), and r > —1: 

r 7^ r = (the limit as r — >■ 0) 



probability mass function = -^^ — ''^-r_i Qk 



k 



k log a 



l-(l-u{l-a)y ^ /X _ log(l-n(l-a)) 



probability generating function gq{u) = jg^ (i_u(i-a)y 9q{u 



log Q 



(where |m| < 1/(1 — a) for the probability generating functions). In our examples we have r = 
1 (geometric distribution) for the Polya Aeppli process, r = (extended truncated negative 
binomial distribution) for the Poisson Inverse Gaussian process, and r = (logarithmic distribution) 
for the Negative Binomial process. 

Example 1 (Polya Aeppli process). We choose 



a = 1 — p and r = 1 for some p S (0, 1) and A > 0. 



The probability mass functions for the two fractional Polya Aeppli processes can be obtained from 
di), (©, (HD and (P with ql"" = (^ij) (1 - p) (see e.g. Example 3.19 in JI^, where P in that 
reference stands for another interesting reference on this process is f2W). 
Alternative formulas for the mixed representation can be obtained from (jl3p and (jl4p by considering 
S{-) as a compound Poisson process with exponentially distributed summands; more precisely, for 
IJL = \/p and /3 = (1 — p)/p, we mean S{-) = Ylk=i'^ ^k^'^ ' ^^s'^e {Y^^^ : n > 1} are i.i.d. random 
variables with density f{x) = /3e~^^l(o.oo)(^)> independent of the (non-fractional) Poisson process 
N^{-) with intensity /x (see e.g. Subsection 11.1.2 in f^], which concerns a slightly more general 
situation with the generalized Polya Aeppli distribution and the Tweedie distribution; in particular 



A(l-P) 

Then we can prove the following result for Example [1] as a consequence of Proposition 13. 1[ 



we recover our case with the Tweedie distribution with parameters (—1, , —^))- 



Proposition 4.1. Let us consider the probability mass functions in ([8]) and ([9]), where M{-) as in 
ExampleUl If f ^ (0, 1), we have 



§tpl{t) - (1 - a)^p^_i(t) = -Xplit) + Xpl_,{t) for all integer k > 1 



^^pUt) = -Xp'^it) 



with the initial conditions Pq{0) = 1 and p^(0) = for all integer k > 1. 
If 1^ € (l,oo), we have 

I - (1 - c^)'^Pk-iit) = Wk{t) - ApLi(i) for all integer k>l 

with the initial conditions Po(0) = 1 and p^(0) = for all integer k>l. 

Remark. The equations in (|18p and (|19p coincide for v = 1. They are equations for the probability 
mass functions of the non-fractional Polya Aeppli process, and we do not know any references where 
they appear. 
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Proof. Firstly, for A; = 0, we have the same equations in Proposition 13. H thus we restrict the 
attention on the case k > 1. If u £ (0, 1), by (jlOp and by taking into account that = {1 — a)^~^a 
(for ah integer h > 1), we have 

%^(t) - (1 - a)^pLi(i) = - ^Plit) + A J](l - af-'apUit) 



df" ' ' ' dt 

1=1 



k-\ 



(1 - a) -\pl_,{t) + A^(l - ar^apl_,_,{t) 

\ i=l / 

k 

XpUt) + Xapl_,{t) + A ^(1 - ay-'apl_,{t) 



1=2 

fc-1 



+ A(l - aM-i(i) - A 5^(1 - af apl_^_i{t) 



i=l 

= -\pl{t) + \pl_,{t). 

If G (1, cxo), the proof follows the same lines; we have to consider (jlip instead of (jlOp . and there 
are suitable changes of signs. □ 

Example 2 (Poisson Inverse Gaussian process). We have 



a = l- r = -\ andX = := |((1 + 2/3)^/2 - 1) for some /3,^ > 0. 



We are not aware of any reference with a formula for the convolution densities g^" in ([7])- So 
([8|) and ([9]) do not give completely explicit formulas for the probability mass functions of the two 
fractional Poisson Inverse Gaussian processes. 

The formulas (jl3p and (jl4p for the mixed representation can be obtained by considering S{-) as an 
Inverse Gaussian process; more precisely S{t) = Y^^0{t) should have density 

Moreover (see e.g. eq. (3.39) in 112]) we have 



E 



exp (-ft ((1- 2/30)1/2 _i)j if0<^ 
oo if 61 > 1 



(21) 

2/3 

We remark that the processes S^{-) (if v E (0,1)) and 5'^(-) (if u € (1,cxd)) for Example [2] can 
be seen as a fractional Inverse Gaussian process. A different fractional version of this process has 
been defined in [15] by subordinating a fractional Brownian motion to an Inverse Gaussian process 
(in analogy with the so-called fractional Laplace motion; see [13]). Some properties of the process 
defined in [TS] are illustrated in [17j. 

Example 3 (Negative Binomial process). We choose 



a = p, r = and A = — logp for some p € (0, 1). 



The probability mass functions for the two fractional Negative Binomial processes can be obtained 
from dSD, ®, (H and §1 "'^^'^ C = pi^^^^^T^' where {\s{k,n)\ : A: € {0, 1, . . . , n}} are 
the unsigned Stirling numbers of the first kind (see e.g. Theorem 6 in \27^ ). 

Alternative formulas for the mixed representation can be obtained from ([T3|) and by considering 
S{-) as a Gamma process; more precisely S{t) should have density 

/(a;) = ^^^ e i-p l(o,oo)(2;). 



9 



5 Processes with two fractional parameters 



In this section we generahze the fractional versions of the processes in Examples [T] and [2] by con- 
sidering a further fractional parameter 77 G (0,1)- Throughout this section we use the notation 
n['^\-) for the process nI^'^{-) in ([5]); then, for each example, we consider independent random 
time-changes of n[^\-) in terms of the processes 5'^(-) in (jl3p and S'^{-) in ()14p . i.e. the fractional 
versions of the subordinator S{-) which appears in the mixed representation of M'^(-) presented in 
Section [3] above. Thus we have 

pl''(t) ■= P{M'''''{t) = k), where M^^'^i-) = N^'^\S'' {■)) , if v e (0, 1); (22) 

pl^^^t) := PiM't^'^it) = k), where M^''^(-) = ivj''^(S'^(-)), ii u e (l,oo). (23) 

We remark that we recover the cases presented above by setting ij = 1. We shall see that the 
governing equations for these probability mass functions have not only a fractional time derivative 
but also the fractional difference operator (1 — B)^ as ([6]). 



5.1 Generalized fractional Polya Aeppli process 

We start with the result. 

Proposition 5.1. Let S{-) = Ef=S'^ ^i'^^ be as in ExampleUl If v ^ (Oil); ^he probability mass 
function in (j22p satisfies the equations 



d^n n.u,,^ 1 / . d-'r 



with the initial conditions Pq'^{^) = 1 and p2''^{0) = for all integer k >1. 
If y ^ (1,00), the probability mass function in (j23p satisfies the equations 



(24) 



dK, 



dt 



(25) 



■^prit) = ^i^-'-i^)pnt) 

with the initial conditions Pq''^{0) = 1 and p^''^{0) = for all integer k >1. 
Remark. The equations in ()24p and (|25p coincide for v = 1. 

Proof. In both cases {v G (0, 1) and G (1, 00)) the initial conditions trivially hold. 
Case G (0, 1). We have 

PiS^it) G G) = PiN^iL^it)) = 0)1g(0) + / fu{x,t)dx (for all Borel subsets G of [0,oo)), 

Jg 

where fi, is the density of the absolutely continuous part of the random variable S'^{t); thus, since 
PiNi^iL^it)) = 0) = E^^ii-fif"), we obtain 

^^(t) = E,^,{-f,tniik=o} + r PiN^^H^'iz)) = k)fAz,t)dz. 

Jo 

Thus we get 

^pTit) = [-i.E^,,{-i.tn + -p{i^+^) mt)) i{fc=o} 



1 / d^ \ r'^ d 
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by (2.4.58) in (2.5) in [5] and an integration by parts (where we take into account that 

P{n['^\l''{0)) = k) = l{fc=o})- The last equation reduces to 

because 

- fiE,^ii-fin + ^ (/" + ^) MO,t) = 0. (26) 



We remark that eq. (|26p can be checked by inspection after noting that: 
by eq. (2.4) in [5\ and some computations; 

by taking into account that we have ^f^^" = ^r^^^f"^ by (2.1.17) in 

Finally we can conclude the proof by using ([6]) with A = 1, and with some computations (where 
we distinguish the cases k = and A: > 1). 
Case u € (l,oo). We have 

P{S''{t) G G) = PiN^iA^/^it)) = 0)1g(0) + / fuix,t)dx (for aU Borel subsets G of [0, oo)), 

Jg 

where is the density of the absolutely continuous part of the random variable S'^{t); thus, since 
PiNf^iA^/^it)) = 0) = e-/^'^'*, we obtain 

?l'^(t) = e-^''^*l{,=o} + / P{n['^\a'/''{z)) = k)Uz,t)dz. 



Pk 

We follow the same lines of the proof for the previous case and we give only some details. Firstly 
we get 



df- 



by (2.2.15) in [11], (5.8) in [5j and an integration by parts (where we take into account that 
P{n['^\a^/''{0)) = k) = l{k=o})- The last equation reduces to 

because 
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We remark that eq. (j27p can be checked by inspection after noting that: /j/(0,t) = * e 
by eq. (5.7) in [5]; 

'!kLf(c^f^ r,,-^V.n _ ( d 

V d{ii^/-) dt- V d{ii^/-) l^^ ^ ^ / 



by (2.2.15) in [TT] (for the equahty = {y.^^'^Y e'^^^"^) and some computations. 

Finally we conclude by using ^ with A = 1 and with some computations (where we distinguish 

the cases /c = and k>l). □ 

We already know that, for r/ = 1, M^''^(-) and M'^'^{-) coincide with the fractional versions in 
(FVl) and (FV2) when we deal with the Polya Aeppli process M(-) as in Example [TJ Then one 
can check that the results in Proposition 14. II can be recovered by considering Proposition 15. II with 
r] = 1. Actually (1 — i3)^p^''^(t) = p^'^{t) —p]j'^i{t) and, if we rearrange the terms in a suitable way, 
the equations (fM|) and ([25]) with r/ = 1 meet (fT8|) and (fT9|) . respectively, by setting fi = ^ = 
and (3 = = as in Example [TJ 

5.2 Generalized fractional Poisson Inverse Gaussian process 

From now on we consider the usual symbol 5{x) for the Dirac delta. Furthermore we consider the 
process S{-) = Y^^p{-) as in Example [21 and we introduce the following notation. 

fy{x,t) = fs'^(t){x) is the density of the random variable S'^{t), for u G (0, 1); (28) 

/i/(x,t) = fs''{t){x) is the density of the random variable S'^{t), for v £ (l,oo). (29) 

The main result in this subsection is Proposition 15.31 which provides governing equations for the 
probability mass functions (p],''^ {t))k>o and (p^''^ (t)) k>o when S{-) is as in Example [2j In its proof 
we refer to Proposition 15. 2^ i.e. a preliminary result for /jy and f,^; thus, in some sense, Proposition 
15.21 concerns the fractional versions of the Poisson Inverse Gaussian process in Example [2l and not 
its generalized version considered here. 

Proposition 5.2. Let and fy he the functions in (j28p and (j29p . respectively. If v £ (0)1); 
have 

/^(^, t) - 2^^ f,{x,t) = 2^^ f,{x, t) (30) 

with the initial conditions fu{x,0) = 6{x) and fuiO,t) = 0. 
If y £ (l,oo), we have 

with the initial conditions /j^(x,0) = 5{x) and fu{0,t) = 0. 

Remark. The equations in (j30p and (j3ip coincide for v = \, and coincide with the equation 
Theorem 3.1 in ]15^ (with the same initial conditions). 
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Proof. We introduce the following notation: g{6), for 6 > 0, is the Laplace transform of g{x); 
g*{(y), for o" > 0, is the Laplace transform of g{t) (we need the second one only for u € (0,1)). 
Furthermore we recall that, by (j2ip . we have 



Ks{-e) = logE[e-^^(i)] = ((1 + 2/30)V2 _ i) . 
Case u G (0,1). We have fu{e,t) = Ey^i{Ks{-ey) by ([H]), and 



^h{e.t) = Ks{-e)~ue,t) (32) 

by eq. (2.4.58) in [TT]. Moreover we have 

^^gj {a)=a\f{cj)-<y^-^g{^) (33) 

by eq. (2.4.63) in [11]. We take the Laplace transforms with respect to t in ([32]) and, by taking 
into account ([33]) with g{-) = fu{0, •) for the left hand side, we get 

then we have 

= (a- - Ks{-e))ft{e,a) - a^-^ = (^a'^ - I + ^(1 + 2/30) _ 

and, after some computations, we obtain 

{a'^ - 2^) - 2^9^ fHO, a) - a^-' (a^ - 2^ + Ks{-e)^ = 0. (34) 

We complete the proof for the case v € (0, 1) by checking that (I30p yields (I34p . This will be done 
in what follows. We take the Laplace transforms with respect to x in (j30p and, after an integration 
by parts where we take into account that fu{0,t) = 0, we have 

then we take the Laplace transforms with respect to t and, by taking into account (j33p with 
g{-) = ■^fy{6, ■) and with g{-) = f,^{6, •) together with some computations, we get 



t=o P 



finally we meet ([M]) because we have g^fuiO, t) = ks{—9) (by ([32]) with t = 0) and fu{9, 0) = 1. 

Case V G (1, oo). We have ]y(d, t) = exp(-(-K5(-0)) V^'t) by and ^fu{0, t) = -Ks{-e)U{e, t) 
and '^^^^fu{0,t) = {Ks{—6))'^fy{6,t) by eq. (2.2.15) in [llj. Then one can obtain the equality 

(^i^/.)(M) + 2|(|^/„)«,,,).,^,/„„.„ (35) 



13 



with some computations. We complete the proof showing that ([3T]) yields ([35]) . We take the Laplace 
transforms with respect to x in (|3ip . and we have 



then we meet (jSSp by considering an integration by parts for the right hand side, where we take 
into account the equality fu{0,t) = 0. □ 

Now we are ready to prove Proposition 15.31 

Proposition 5.3. Let (p^''^(i))fc>o (^i^d {p],''^ it))k>o be the probability mass functions in (j22p and 
(j23|) . respectively, with S{-) = ^^,/3(-) as in Example\^ If v ^ (0)1)? have 

^§;Pnt) - 20frW = 2^(1 - BTpnt) for aU ^nteger k > (36) 

with the initial conditions Pq''^{0) = 1 and p^''^(0) = for all integer k >1. 
If € (l,oo), we have 

^^Pr(0 + 2^^PrW = 2y (1 - BrPVit) for aU ^nteger k > (37) 

with the initial conditions Pq''^{0) = 1 and pl''^{0) = for all integer k > 1. 
Remark. The equations in (|36p and (j37p coincide for u = \. 



Proof. In both cases {v G (0, 1) and v € (l,oo)) the initial conditions trivially hold. 

We start with the case v G (0, 1) and we consider the function fi, in (i28]l . Firstly we have 

p^'^(t) = f P{n['^\z) = k)f^{z,t)dz (for aU integer k > 0). 
Jo 

Then we get 

^^pYH) - 2^^pY{t) = / P(ivf')(z) = k) ( ^^f,(z,t) - 2^?^f,{z,t)] dz 
df dt" pdt^^'' Jo ^ 1 ^ ^ ^dt'^dt''-'^ ' ^ pdt" ^ ') 



P{N[''\z) = k)2^—f,{z,t)dz 



by ([30P in Proposition 15. 2| then 



=2^i-{l - B)' r P{Nf{z) = k)Uz,t)dz 

P Jo 



by an integration by parts and by taking into account fu{0,t) = (for the first equality) and by 
([6]) with A = 1 (for the second equality) . Thus (j36p holds and this completes the proof for the case 
u € (0, 1). The proof for the case u € (l,oo) follows the same lines; we get ([57P by considering a 
suitable change of sign and, in particular, ([3T]) instead of ([30]) . □ 

We already know that, for r] = 1, M^''^[-) and M^''^{-) coincide with the fractional versions in 
(FVl) and (FV2) when we deal with the Poisson Inverse Gaussian process M(-) as in Example 
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[H Then one can check that, for k = 0, the equations in Proposition 15.31 can be obtained from the 
equations in Proposition 13.11 adapted to this process. Actuahy, for u £ (0, 1), we have 



^pl^^it) = -\p,,pYit) = -^((1 + 2/3) V2 _ i)pY{t) 



by (jlOp in Proposition 13.11 then we obtain 



with some computations, and this meets ([36]) in Proposition 15.31 because (1 — BYp^j^'^{t) = p]^'^{t) 
Pk-i{t) and pi'i(t) = for all t>Q. For v G (l,oo) we have similar computations with suitab 
changes of signs (we have pT]) and ([37|) in place of ([10]) and ([36]) . respectively). 
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